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Abstract. Cosmic ray (CR) air showers, detected via the air fluorescence technique, 
are reconstructed in part using functions that parameterize the longitudinal profile of 
each shower. The profile parameterization yields the position of shower maximum, 
Xmax hi gm/cm^, which is sensitive to the incident CR particle type: e.g. p, C/N/0, 
Fe or 7. The integral of the profile is directly related to the shower energy. The Pierre 
Auger fluorescence reconstruction uses the Gaisser-Hillas (GH) 4-parameter form[5]. 
The HiRes group has used both the Gaisser-Hillas form and a 3-parameter Gaussian 
in Age form[7]. Historically analytic shower theory suggested yet other forms; the best 
known is a 3-parametcr form popularized by Greisen[6]. Our work now uses the shower 
full width half-maximum, fwhm, and shower asymmetry parameter, /, to unify the 
parameterization of all three profile functions. Furthermore shower profiles expressed 
in terms of the new parameters: fwhm, f have correlations greatly reduced over e.g. 
Gaisser-Hillas parameters Xq,X. This allows shower profile reconstructions to add 
constraints (if needed) on the mostly uncorrelated parameters fwhm, f. 
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1. Introduction 

Following the first observation of extensive air showers by the air fluorescence 
technique[l], the Fly's Eye experiment pioneered analysis techniques[2, 3, 4] that 
are largely unchanged to this day. Most relevant to this paper are the analysis 
techniques needed to exploit the feature of fluorescence detectors to observe directly 
the development profile of air showers in the atmosphere. The shower profile provides 
a direct measurement of the depth of shower maximum, Xmax in gm/cm^, and provides 
an almost model- independent measurement of the cosmic ray energy [4]. 

Historically, the shower profile was parameterized [2] using three parameter 
Gaussian or four parameter Gaisser-Hillas[5] (GH) functions to characterize the shower 
brightness versus shower depth, X(gm/cm^). The Gaussian parameterization was 
appealing, as it depended on distance from shower maximum, X — Xmax, and depended 
on one parameter to characterize the shower width. As the precision of fluorescence 
detector measurements improved the Gaussian parameterization was discarded, as it 
did not model the true shape of shower profiles that are asymmetric in shower depth 
with respect to shower maximum. The use of the four parameter GH function was often 
limited by the realities of data signal uncertainties or by the limited fraction of the 
shower profile observed for any given shower. 

More recently the HiRes/MIA collaboration [7] analyzed a composite^ shower profile 
using three functional forms: Gaisser-Hillas[5] (noted previously), a form motivated 
by analytic shower theory and popularized by Greisen[6] (Greisen), and a new form 
motivated by the observation that showers were rather symmetric when plotted in 
shower age, termed Gaussian in Age (GIA)[7]. Interestingly, the three functions were 
almost indistinguishable and described the data comparably well. The paper found a 
strong correlation between two of the GH parameters and suggested that a simpler (3 
parameter) form should be adequate to describe the shower profiles. A closely related 
comparison of GH and GIA functions to Monte Carlo simulated showers [8] generated 
using C0RSIKA[9], again found comparably good agreement using GH or GIA profile 
functions. Interestingly the Monte Carlo study [8] observed the near equality of the 
width at half-maximum, fwhm, of proton, iron and photon showers but did not exploit 
this fact. 

Motivated by the history above, in this paper we write the GH, GIA and Greisen 
functions in terms of 4 parameters: the intensity at shower maximum, Nmax, the depth 
of shower maximum, X^ax, and two other parameters that in the case of GH are the 
parameters Xq and A§. In all cases the profile functions can then be re-expressed in 
terms of a physical distance A = X — Xmax, a composite parameter, W = Xmax — Xq, 

t Showers observed by the HiRes prototype/Ml A detectors were ahgned to have a common shower 
maximum. The showers were then averaged to obtain a composite shower profile that was then used 
in their analysis. 

§ While the conventional versions of GIA and Greisen are recovered setting Xq ~ 0, the introduction 
of Xq allows for a more symmetrical comparison to GH as well as allowing these forms to be applied 
to e.g. neutrino induced showers. 
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and a (forth) dimensionless parameter, ^. We will now show that W depends both on 
shower width, e.g. fwhm, and on a shower asymmetry parameter, /. The parameter, 
^ depends only on the shower asymmetry, /. 

2. Profile functions 

We now summarize the three profile functions: GH, GIA and Greisen. 

2.1. Gaisser-Hillas 

The Gaisser-Hillas function [5] is conventionally written: 

N{X)gH = N,nax (y ' ^ ~ 

If we introduce the new variable: W = X^ax — and the distance from shower 
maximum: A = X — Xmax, the GH function can then be written as: 

N{X)gh = N^ax (1 + ^)'^ e-x (1) 

It is useful to introduce two dimensionless quantities: e = and ^ = Then a 
minimal form, useful for comparisons, is obtained: 

N{X; N^ax, Xmax. ^0, \)gh = N^^x (1 + ef e-(^ (2) 

This form emphasizes the role of a physical quantity, the distance from shower maximum, 
A, in comparison to e.g. an unphysical quantity, X — Xq. Furthermore A is scaled by 
a parameter W potentially resolving a tension between X^ax and Xq in a parameter 
optimization to fit experimental shower profiles. 

While the new variables, ^, are less intuitive than the standard GH parameters 
Xq, a, we will show that they are related to two very intuitive parameters: the shower 
width, fwhm and a dimensionless shower asymmetry parameter, /. 

For later reference the integral of the GH function[18] is given by: 

/"OO 

/ N{X)GHdX = N^ax W r^^+'^ e« + 1) (3) 
Jo 

where F is the Euler Gamma function. 

2.2. Gaussian in Age 

The Gaussian in Age function[7] is conventionally written: 

N{X)gia = Nmax e-^ (^)' 



where "s" is the shower age and s = 1 at shower maximum. Motivated by the definition 
of the GH function, the conventional definition for Sconv = vJ^v — is then modified by 
replacing: X X — Xq and Xmax Xmax—X^. Then with the same definitions as used 
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for the GH function (above), the redefined age variable is: s(e) = j^^- Note that this 
definition preserves s = 1 at shower maximum and s = 3 in the hmit of X >> X^ax, but 
differs from "Scon?;" at other shower depths. As with the GH function, a minimal form 
for the GIA function is obtained: 

N{X;Nmax,X,nax,Xo,a)GiA = N^ax e'^^^^' (4) 

2.3. Greisen 

The Greisen function[6] is conventionally written[6, 7]: 

]\T(f\^ . — M pi t(l-|'™(s)) - tmax ) 

jGreisen -'■"max ^ 

where t is the depth in the shower in radiation lengths, tmax is the depth of shower 
maximum and "s" is shower age. Following the modification of GIA: t — t — to and 
tmax — tmax — to, then Converting from depth in radiation lengths to depth in gm/cm^, 
t = 3g~7 ' and using the same notation as for GH and GIA (above) we obtain the minimal 
form: 

N(Y- N X X r> \ - N e(l-fin(s(e))) - |«n(s(e)) ) ^ 

ly l^max-, ^max-, ^0^P-iG.7jGreisen — ^^max ^ ^ ^36. 7 

where the conversion constant "36.7gm/cm^" (per radiation length) is treated as a 
parameter [7], Pae.y li- 
lt is interesting to note that Eqn. 2, 4 and 5 depend only on e = and one 
additional dimensionless parameter: ^ = a or As noted above, we will now 
show that W depends on both shower width and asymmetry: fwhm, f and ^ depends 
only on shower asymmetry: /. 



2.4- Profile parameters: fwhm and f 

While there are now some similarities between the three profile functions, a simple 
connection between the profile parameters and the appearance of a shower does not 
exist. A natural choice is to use the width of the shower (at e.g. half-height), fwhm 
and an asymmetry parameter, /. Using the profile in Fig. l(left), fwhm = C + TZ 
and / = C/{C + 7V). Thus fwhm has dimensions (gm/cm^) and / is dimensionless. 
Furthermore, because real and simulated shower profiles rise more steeply than they fall 
/ will always be less than 0.5. 

CORSIKA simulations [8] had suggested that shower fwhm values were rather 
similar for p, iron and 7 showers. Results of CONEX[10, 11, 12] simulated showers 
are shown in Fig. 2 for fwhm and for the asymmetry parameter, /, for shower energies 
relevant to the Auger [13], HiRes[14] and TA[15] experiments. These simulations are 
consistent with the CORSIKA results and suggest that for energies ~ 10^^'^eV and for 
all primaries the fwhm is in the rather limited range: 475 gm/cm^ ^ fwhm ^ 575 
gm/cm^ and asymmetry in the range: 0.44 ^ / ^ 0.47. Furthermore the CONEX 

II The original Greisen functional form is simply regained by setting: Xq ~ and P36.7 = 36.7. 
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Model shower profile: Xmax = vao gra/cm^, fwhra = 500 gm/cm', f=o.45 Variations of Gaisser-Hillas Profile with asymmetry (f) 




Figure 1. Left: Plot of a representative Gaisser Hillas (GH) shower profile defined 
by: Xmax — 700 gm/cm^, fwhm = C + TZ = 500 gm/cm^ and shower asymmetry 
f = £/{C + TZ) = 0.45. The vertical axis is the shower intensity {Nf.) normalized to 
the intensity at shower maximum (Nmax) and the horizontal axis is the shower depth 
{X) in the atmosphere in gm/cm^. Right: GH shower profiles with X,nax = 725 
gm/cm^, fwhm ~ 525 gm/cm^ and three different values of asymmetry: / — 0.44, 
0.45 and 0.46. The corresponding GH parameters are: {Xq ~ 2.2 gm/cm^, A ~ 68.3 
gm/cm^), {Xq = —145 gm/cm^, A = 56.9 gm/cm^) and (X^ — —365 gm/cm^, A = 45.5 
gm/cm^) respectively. 

simulations suggest that the asymmetry parameter, /, may provide some discrimination 
in primary composition [16, 17], i.e. between different types of primary particles. Finally 
both fwhm and asymmetry, /, show a slow, approximately logarithmic increase with 
increasing shower energy. 



Stody Conex Stiowers: profile fwtim Study Conex Stiowers: profile asymmetry, f 




16.5 17 17.5 18 18.5 19 19.5 20 20.5 ' 16.5 17 17.5 18 18.5 19 19.5 20 20.5 

log,o shower energy log^g shower energy 



Figure 2. The results of CONEX[10] simulations of proton, iron and 7 initiated air 
showers are shown in blue, red and green respectively. Left: Plot showing the shower 
widths (fwhm) on the vertical axis versus the shower energies on the horizontal axis. 
The error bars give the RMS of the simulated showers in each energy bin. Right: Plot 
showing the shower asymmetries (/) on the vertical axis versus the shower energies on 
the horizontal axis. 
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2.5. Relating old and new shower parameters 

To relate the old to the new shower parameters, the shower profile is evaluated at a 
fraction h of maximum. N{X) = h Nmax- To illustrate the procedure we choose the GH 
profile as written in Eqn. 1. Then solving for the values of A that satisfy A^(^) = h N^ax 
one obtains: 

W A A 

This relation is true for two values of A, see Fig. l(left): = {—f{h)) Width{h) and 
Aji = (1 — /(/i)) Width{h) where jA^I + A-ji = Width{h) is the width of the shower 
profile curve at fractional height h^. Following the familiar convention to describe a 
profile at half-height, we define fwhm = Width{h = 0.5) and / = f{h = 0.5). 



IS GH, Gaussian in Age GIA, and Greisen Profile 
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Figure 3. Left: The curves show the ratio function. plotted on the vertical 

axis versus the shower asymmetry, /, on the horizontal axis. As discussed in the text 
the ratio function: R = = x^^^™Xo depends only on the shower asymmetry. 

The ratio functions for the different shower profiles arc plotted as: GH (solid). GIA 
(dashes) and Greisen (dots). Right: The asymmetry /-dependent terms in Eqn. 7, 

^ hOT^^^^^' plotted on the vertical axis versus the shower asymmetry, /, on 

the horizontal axis. 



The equations for A^ and A-ji can be combined to obtain one (GH profile specific) 
relationship: 

l-fR 



ln{ 



l + il-f)R' 



R = 



(6) 



where R = depends only on the shower asymmetry, /. Eqn. 6 is solved numerically 
to obtain the ratio function: R = R{f)', the result is shown in Fig. 3(left). 

The solution to Eqn. 6 relates W to the shower fwhm: W = Finally with 

the ratio function R{f) and W known, the GH parameters Xq and A are given by: 



^ While both Width{h) and f{h) depend on the choice for h, the shower shape and the parameters 
describing the shower shape, e.g. GH parameters Xq, A, do not. 
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X 



— Xmax — W 

W [ ln{l - fR) + fR 



ln{h = 0.5) 



It is important to note that ^ = x depends only on the asymmetry parameter, /. 
Thus the calorimetric shower energy, given by: 

POO 

EThZer = / N{X)gh dE/dx dX 
Jo 

is to a good approximation[19], see Eqn. 3: 

EThoZer = < dE/dx > N^a, fwhm C- R(ff^^ (7) 

where < dE/dx > ~ 2.19 MeV/particle/(gm/cm^) [19] and N^ax is the number 
of charged particles at shower maximum. The dependence on shower asymmetry, / 
is shown in Fig. 3 (right) and is ~ 1% over the range of asymmetries predicted using 
CONEX, see Fig. 2 (right). Thus to a good approximation the shower calorimetric energy 
is proportional to: Nmax fwhm; see Fig. 4. This relation emphasizes the importance 
of the shower width, fwhm, as errors in determining the width directly propagate into 
errors in estimating the shower energy. This relation also shows that mis-measuring the 
asymmetry has almost no effect on the shower energy. Eqn. 7 and Fig. 3 (right) do imply 
a small correlation coming from the calorimetric shower energy constraint: namely that 
more asymmetric showers (smaller values of /) should have slightly smaller values of 
fwhm. 

To relate conventional GIA parameters, Xq, a, to new parameters, fwhm, f, we 
follow the same procedure that was employed for the GH function. For the GIA profile, 
the simple mathematical form of the GIA function (see Eqn. 4) results in an analytic 
(GIA profile specific) solution for the ratio function R{f)'. 

R(f) - 

which is plotted in Fig. 3(left). As with GH, W = The conventional GIA 

parameters are then given by: 

fR 



a 



{3-fR) ./zM!}=m 



2 

Finally to relate conventional Greisen parameters, Xo,p36.7, to new parameters, 
fwhm, f, the first step is to solve for the values of A that satisfy N{X) = h N„ 



''max- 
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Figure 4. Scatter plot showing CONEX simulations of proton, iron and 7 initiated air 
showers shown in blue, red and green respectively. The points report the shower fwhm 
in (gm/cm^) on the vertical axis versus the ratio of shower quantities: N max I E shower 
in (particles/GeV) on the horizontal axis. See Eqn. 7 for motivation. All showers 
had an energy near 10^^'^eV. The relative (horizontal) displacement of the iron 
: p : 7 showers is related to increasing missing energy (in muons and neutrinos) in 
iron versus p versus 7 showers. The points with error-bars reflect the average of the 
proton points in the scatter plot. 



As with the GH function, this relation is true for two values of A: and At^. The 
two resulting equations can be combined into one relationship that relates R = 
to the shower asymmetry, /. The solution, found numerically, is plotted in Fig. 3(left). 
Once the ratio function -R(/) is determined, then the Greisen parameters, Xo,p36.7 are 
given by: 

-^0 = ^max — W 

W 3 3 

= ln{h = 0.5) f-^"^^' - 2^"^'^^) - 2^"^'^)] 

where sc = ^i^) is the (redefined) shower age (see Sect. 2.2 above). 
3. Profile comparisons 



Now that the GH, GIA and Greisen functions can be expressed in common parameters: 
fwhm, f, it is easy and instructive to compare their shower profiles. An example 
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Gaisser Hillas, Gaussian in Age, and Greisen Profiles: Xmax = 725 gm/cm^, fwhm= 525 gm/cm^ Gaisser Hillas, Gaussian in Age, and Greisen Profiles: Xmax = 725 gm/cm^, fwhm= 525 gm/cm^ 




Figure 5. Left: Plot of shower profiles for GH (solid), GIA (dashes) and Greisen 
(dots) profile functions for typical shower profile parameters: X^ax = 725 gm/cm^, 
fwhm = 525 gm/cni^, and / = 0.45. The vertical axis is the shower intensity 
normalized to the intensity at shower maximum and the horizontal axis is the shower 
depth in the atmosphere in gm/cm^. Right: Same shower profiles now with a semi-log 
plot to show the small differences between the functions. 



comparison, shown in Fig. 5, finds that GH and Greisen functions are almost 
indistinguishable while (both) being very similar to, but systematically below the GIA 
profile for shower depths well away from shower maximum. This result is consistent with 
HiRes/MIA[7] and HiRes[8] observations that the profiles describe experimental (and 
simulated) shower data comparably well near shower maximum. Quantitatively, the 
small but systematic profile function differences result in shower calorimetric energies 
evaluated using the GIA function being ~ 1% larger than those evaluated using GH or 
Greisen forms. 

As noted in Sect 2.4, the shower profiles away from shower maximum may provide 
some discrimination on cosmic ray composition [16, 17] (cf. Fig. 2(right)) and on whether 
the GIA [20] versus GH (or Greisen) profiles provide somewhat better parameterizations 
of actual shower profiles. While interesting, this is not the thrust of this paper and will 
not be discussed further. 

A remaining issue is how best to apply four parameter shower profile functions to 
less than ideal data. In these circumstances the tendency is to fix, or constrain, one or 
more of the parameters. An equivalent point of view is the assumption that Xq = 
gm/cm^ in the GIA or Greisen function and/or that the parameter p^Qj = 36.7 gm/cm^ 
(per radiation length) in the Greisen function. 

To study this issue, shower profile parameters from CONEX showers are shown in 
Fig. 6 for proton showers with energies ~ 10^^'^eV. For each simulated shower, the shower 
fwhm and asymmetry, /, are determined. The values are plotted in Fig. 6(lower right). 
Following Sect. 2.5 the fwhm and / values are used to determine the corresponding 
parameters for GH, GIA and Greisen profiles. The results are plotted in the other 
quadrants of Fig. 6. The conventional shower shape parameters: (Xo,A), {Xo,a), 
and (^0)P36.7) show significant correlations while {fwhm, f) do not. In particular 
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Figure 6. Scatter plots of Top Left: GH parameters (Xq, A), Top Right: GIA 

parameters (Xq, cr), Bottom Left: Greisen parameters (^o: P36.7) foi' proton and 
Bottom Right: shower {fwhm, /) for proton (blue), iron (red) and 7 (green) showers 
near 10^^'^eV simulated with CONEX. hi all cases the first parameter is plotted on the 
vertical axis and the second parameter on the horizontal axis. The points with error- 
bars reflect the average of the (proton) points in the scatter plot. The plot statistics 
are: 411 proton, 440 iron and 462 photon CONEX showers. 



the correlation coefficients for GH, GIA and Greisen parameters are in the range: 
0.76 ~ 0.91 for CONEX[10, 11, 12] proton and iron showers while the correlation 
coefficients for {fwhm, f) are ~ 0.12 for proton showers and ~ —0.11 for iron showers. 

While the GH correlation results appear consistent with the observation of 
parameter correlations in Ref.[7], our analysis does not (then) support the use of 
GIA (or Greisen) shower proffies with the Xq parameter set to zero[7]. Rather than 
the correlation being evidence that three parameters are sufficient to describe shower 
profiles, it is more likely that GH, GIA and Greisen parameters need to fulffil the 
calorimetric energy constraint of Eqn. 7"*". As noted above, this constraint is almost 
insignificant with the choice of variables {fwhm, /). In summary the {fwhm, f) shower 
shape parameters appear less correlated that the other shower parameter pairs (while 
still assuming a wide range of values) . 

That said, it is likely that shower details lead to correlations with any choice of 

+ An almost equivalent statement, is that e.g. {Xq,\) must be compatible with the full with half 
maximum of the shower. 
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shower parameters. In that regard, this study is compatible with shower reconstructions 
using partial constraints [18]; we recommend that the constraints be applied to fwhm 
and / rather than the GH parameters Xq and A. 

4. Conclusion 

This analysis has studied three functions used to parameterize the shower profiles of 
ultra-high energy cosmic rays: Gaisser-Hillas (GH), Gaussian in Age (GIA) and an 
analytical form popularized by Greisen (Greisen). As part of this analysis, we have 
generalized the GIA and Greisen functions to 4-parameters consistent with the GH 
function. To then study the relation between the three shower parameterizations, and 
to use more physically descriptive parameters, this work uses the shower full width half- 
maximum, fwhm, and shower asymmetry parameter, /, to parameterize the shape of 
the shower profiles. For profiles with the same {fwhm, /), the GH and Greisen shower 
profiles are essentially identical and systemically less than GIA for shower depths away 
from shower maximum. Of the three functions, GH is most convenient as the integral of 
the GH profile is an analytic function. Monte Carlo simulated air showers using CONEX, 
and parameterized in terms of the new parameters: {fwhm,f), have correlations 
(between those parameters) greatly reduced over the standard parameterizations e.g. 
Gaisser-Hillas parameters: {Xq,X). This allows shower profile reconstructions to add 
constraints (if needed) on the mostly uncorrelated parameters fwhm, f. The CONEX 
shower simulations suggest that the shower asymmetry parameter, /, may have some 
sensitivity to the incident cosmic ray particle type: e.g. p, C/N/0, Fe or 7. 
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